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Abstract 

This paper deals with the modeling of the interstellar hydrogen atoms (H atoms) 
distribution in the heliosphere. We study influence of the heliospheric interface, that is 
the region of the interaction between solar wind and local interstellar medium, on the 
distribution of the hydrogen atoms in vicinity of the Sun. The distribution of H atoms 
obtained in the frame of the self-consistent kinetic-gasdynamic model of the helio- 
spheric interface is compared with a simplified model which assumes Maxwellian 
distribution of H atoms at the termination shock and is called often as 'hot' model. 
This comparison shows that the distribution of H atoms is significantly affected by the 
heliospheric interface not only at large heliocentric distances, but also in vicinity of 
the Sun at ~ 1 — 5 AU. Hence, for analysis of experimental data connected with direct 
or undirect measurements of the interstellar atoms one necessarily needs to take into 
account effects of the heliospheric interface. In this paper we propose a new model 
that is relatively simple but takes into account all major effects of the heliospheric 
interface. This model can be applied for analysis of backscattered La-alpha radiation 
data obtained on board of different spacecraft. 
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Introduction 



The Solar system is moving through a partially ionized Local Interstellar Cloud (LIC) 
with relative velocity of Vuc = 26.4 km/s (Moebius et al. 2004). The main component 
of chemical composition of the LIC is atomic hydrogen (about 90%). Interstellar atoms 
of hydrogen (H atoms) can deeply penetrate inside the heliosphere due to their large free 
path. These atoms effectively interact by charge exchange with protons of the interstellar 
cloud and the solar wind. Direct and indirect measurements of H atoms at the Earth orbit 
can provide substantial information on the properties of both the heliospheric interface and 
the Local Interstellar Medium (LISM). Theoretical modeling of distribution of the inter- 
stellar H atoms inside the heliosphere is fairly important to understand physical processes 
in the heliosphere and for correct interpretation of experimental data. At present main 
information on the properties of interstellar H atoms inside the heliosphere are obtained 
from backscattered Ly-alpha radiation measurements obtained by SOHO, HST, Voyager 1, 
Voyager 2, Pioneer- 10, Cassini, Ulysses, Galileo and other spacecraft. The diagnostics of 
the H atom distribution by backscattered Lyman-alpha radiation is possible because spec- 
tral properties of the radiation considerably depend on distribution function of H atoms 
(Quemerais 2006). 

During several last decades the classical hot model is generally employed to analyze 
backscattered Ly-alpha radiation data (Wu and Judge 1979; Lallement et al. 1985). With 
this model one can obtain a solution of the Boltzmann's equation for the velocity distribu- 
tion function of H atoms. However, the model has several crude assumptions that make 
the model less appropriate for correct analyzes of experimental data. These assumptions 
are: 

• Maxwellian distribution function is taken at infinity (i.e. in LISM) as the outer 
boundary condition: 

lim f{wH, r) = /m(wh). 

r-^co 

This condition means that no effects of the heliospheric interface can be taken into 
account in the model. The heliospheric interface effects can be taken only by artifi- 
cial changes of the interstellar bulk velocity and temperature. 

• The model is stationary. This means that no variations of the solar wind parameters 
with the solar cycle can be taken into account in the frame of the model. 

• The model is axisymmetric. This means that the effects connected with the helio- 
latitudinal variations of the solar wind's parameters or interplanetary and interstellar 
magnetic field can be taken into account. 
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Figure 1 : A: Qualitative structure of the heliospheric interface. HP is the tangential dis- 
continuity (heliopause) separating the fully ionized solar wind from the interstellar plasma 
component, BS is the bow shock, and TS is the heliospheric termination shock. The dashed 
lines show the trajectory of a solar wind atom of populations 1 and 2 and of an interstellar 
atom of populations 3 and 4. B: cylindrical coordinate system that is employed in this 
work, z-axis is the axis of symmetry, it is directed opposite to upwind; curve 1 is the outer 
boundary; curve 2 is the heliospheric termination shock. 

In recent modifications of the hot model (Lallement et al. 1991, Bzowski & Rucin- 
ski 1995) the boundary conditions are taken far from the Sun (particularly at sphere with 
radius 80-90 AU), where solar gravitation, radiation pressure and photoionization are neg- 
ligibly small. Parameters (i.e. number density, average velocity and average temperature) 
of the Maxwellian velocity distribution function of H atoms at the boundary were used 
as free parameters and chosen in order to get optimal agreement between theoretical and 
experimental results. Costa et al. (1999) determined the temperature and velocity of inter- 
stellar atoms in the distant heliosphere (at 50 AU) on basis of analysis SWAN/SOHO data. 
It was shown that these values are significantly different from local interstellar parame- 
ters known from the measurements of interstellar helium (Moebius et al. 2004). These 
discrepancies are explained by heating and deceleration of hydrogen atoms's flow during 
their motion through the heliospheric interface (see Baranov and Malama 1993). Fig.[lK 
demonstrates a qualitative picture of the heliospheric interface region. The heliopause is 
a contact discontinuity, which separates the solar wind and the interstellar plasma com- 
ponent. The termination shock (TS) is formed due to the deceleration of the supersonic 
solar wind. The bow shock (BS) may also exist if the interstellar plasma flow is super- 
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sonic. Four regions can be distinguished: the supersonic solar wind (region 1); the solar 
wind flow between the TS and the HP (region 2 or the inner heliosheath); the disturbed 
interstellar plasma component flow (region 3 or the outer heliosheath); the undisturbed 
interstellar gas flow (region 4). 

The first self-consistent kinetic-gasdynamic model of the solar wind interaction with 
the two-component (neutral atoms and plasma) LIC had been developed by Baranov and 
Malama (1993). In this model the kinetic equation for the distribution function of H atoms 
is solved self-consistently with the Euler equations for the charged component. The main 
physical process considered in the model is the charge exchange process of the H atoms 
with protons: + H = H + Atoms newly created by charge exchange have the 
velocity of their ion counterparts in charge exchange collisions. Therefore, the velocity 
distribution of these new atoms depends on the local plasma properties in the place of 
their origin. It is convenient to distinguish four different populations of atoms depending 
on the region in the heliospheric interface where atoms were formed. Population 4 are the 
primary interstellar atoms, which suffer substantial filtration in the heliospheric interface. 
The mean free path of the atoms with small velocities is smaller as compared with the 
fast atoms. Therefore, slow atoms are more processed by charge exchange (and ionization 
processes in general). This kinetic effect called the selection effect, results in asymmetry 
of the distribution function of population 4 at the TS (see Izmodenov et al. 2001). Popu- 
lation 3 of H atoms in the heliosphere are the secondary interstellar atoms created in the 
outer heliosheath. Population 2 are the atoms created in the inner heliosheath between 
the heliopause and the termination shock. The number density of these atoms is signifi- 
cantly smaller as compared with the number densities of primary and secondary interstel- 
lar atoms. The population 2 has no major importance for interpretation of backscattered 
Lyman-alpha radiation experiments, but has major importance for the Interstellar Bound- 
ary Explorer (IBEX) spacecraft measurements that are specifically designed for measure- 
ments of this component. Population 1 consists of the atoms created in the region 1 of the 
supersonic solar wind. It can be noted, atoms of this population have velocity ~ 400 km/s, 
that leads to a large doppler shift in line of Ly-alpha. Hence, atoms of population 1 do not 
add backscattered photons in Lyman-alpha. 

The mean free pathes with respect of charge exchange are comparable or larger than 
the characteristic size of the heliospheric interface for all introduced populations of H 
atoms (see, e.g., Izmodenov et al. 2000). While the interstellar H atoms passing the inter- 
face before entering inside the supersonic solar wind region their distribution function is 
disturbed significantly due to the process of charge exchange. It was shown in Izmodenov 
et al. (2001), that distribution function of all populations is very different from Maxwellian 
at any given point of the heliospheric interface. Actually, the proof of the non-maxwellian 
behavior of the interstellar H atoms in the heliosphere was done by Baranov et al. (1998) 
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who have shown that Tr and Tj- temperatures are different. The temperatures are defined 
as: 

Tr{v) = / /(r,w) ■ iVR{r) -WR^dw 

= / ■ (V,(r) -w.fdw, 

k-nH[r) J 

here /(r, w) is the velocity distribution function, r, w is the radius-vector of atom and its 
velocity, respectively. n//(r) is the number density, Vr and are components of the bulk 
velocity in the cylindrical system of coordinate (with z axis as the axis of symmetry). niH is 
the mass of atomic hydrogen, k is the Boltzmann constant. Significant difference between 
values of component Tr and shown by Baranov et al. (1998) has proven that there is an 
asymmetry of the distribution function in different directions. 

In this work we explore the influence of the heliospheric interface on the distribution of 
the interstellar H atoms inside the heliosphere. To do this we compare the results obtained 
in the frame of the hot model with the results obtained in the frame of the heliospheric 
interface model (Baranov and Malama, 1993). The comparison illuminate the role of the 
heliospheric interface on the distribution of H atom in the supersonic solar wind region. We 
also present an advanced hot model that combines the simplicity of the classical hot model 
with a possibility to take into account the heliospheric interface effects on the distribution 
of the H atoms. 

Model 

We study the velocity distribution function and its momentums of the interstellar H atoms 
inside the heliospheric termination shock (i.e. in the region 1 introduced above). In the 
considered model we obtain the distribution function inside the sphere with center at the 
Sun and radius 90 AU. This sphere serves as the outer boundary. It is located inside the he- 
liospheric termination shock and in upwind direction it is close to the TS (see Fig.[T]B). We 
use the following notations: z-axis is directed in the opposite direction of the interstellar 
flow, 7?-axis (in cylindrical coordinate system) is directed in the direction perpendicular 
to the z-axis. 9 is the angle counted from positive direction of the z-axis. In this pa- 
per we restrict ourselves to the stationary axisymmetric problem (z is a symmetry axis). 
This corresponds to stationary and spherically symmetric flow of the solar wind. How- 
ever, the method and numerical code that was developed by us can also be used for the 
non- stationary and three-dimensional cases without any changes. 
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In the solar system atom is affected by the solar gravitational force Fg and the solar 
radiation pressure F^-^^. These forces act in opposite direction and both are proportional to 
l/r^. We can introduce dimensionless parameter 

which determines a balance between the solar gravitational and radiation pressure forces. 
Resulting force acting on unit of mass is represented as: 

here G is the gravitational constant. My is mass the Sun, r is the radius-vector. In general 
case parameter ji depends on time and radial component of velocity v^. However, for 

the purposes of this paper we assume that ji = 1.258 and this value does not change 
with time or atom velocity. Therefore, we will restrict the study to the case when the result- 
ing force F is repulsive. The chosen value of ju is obtained by averaging of experimental 
data during one solar cycle. The kinetic equation for the velocity distribution function of 
H atoms can be written as: 

df(r,w) <?/(r,w) ^ df(r,w) ^. , (9F or x x 
at or ow ow 

Here r,w are the radius-vector and the velocity of H atom, respectively; /(r,w) is the 
velocity distribution function of H atoms; F is the force as described above. The last term 
in the left part of this equation appears in the 'exotic' case, when the effective force de- 
pends on the radial component of the H atom velocity (see Longmire 1963). For case 
when n does not depend on the velocity of atom this last term is equal to zero. The right 
part of the kinetic equation includes the loss of atoms due to the charge exchange process, 
photoionization (H + hv = H'^ + e) and electron impact ionization (H + e = H'^ + 2e). 
Coefficient /3(r) is the effective ionization rate: j8(r) = ^exir) + Pphif) + ^eiir) , where 
Pex{f),Pph{f) Pei{f) ^rc the ratcs of the charge exchange, photoionization and electron 
impact ionization, respectively. It is assumed in the model that the rates of charge ex- 
change, photoionization and electron impact ionization decrease with distance from the 
Sun as ~ 1 /r^, where r is the heliocentric distance. This assumption deals with the fact 
that the rates Pex, Pei and jSp/, are proportional to the number density of the solar protons, 
electrons and photons, respectively, and in the stationary case these number densities are 
proportional to ~ 1/r^. Thus: 
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Figure 2: Comparison of the results of the hot and Baranov-Malama models. Top raw (A) 
shows parameters for the primary interstellar atoms: Tr and temperatures in the upwind 
and crosswind directions. Bottom raw (B) shows parameters of the secondary interstellar 
atoms: number density (in dimensionless units) in upwind and crosswind, temperature Tr 
in upwind and temperature in crosswind. Number density for primary H atom compo- 
nent is not presented because the curves of number density coincide for different models 
that makes the comparison not informative. 

Solid curves correspond to the Baranov-Malama model; dashed curves correspond to hot 
model with Maxwellian distribution at the outer boundary. 



Here = 1 AU is the distance from the Earth to the Sun. Symbol E shows the values 
at the Earth's orbit. In this paper we assume = 6.2 x 10^^ s^^ This magnitude was 
obtained by averaging of available relevant experimental data over the solar cycle. 

Generally speaking, new hydrogen atoms (population 1) are created inside the com- 
putational domain due to the charge exchange process. These atoms have properties (in 
particulary the velocity and temperature) of the supersonic solar wind. Therefore, a photon 
scattered by such an atom has large doppler shift and does not contribute to the backscat- 
tered Ly-alpha spectrum. Thus, population 1 is not important for the present study and the 
charge exchange process results simply in losses of other populations of H atoms inside 
the helio sphere. 

The kinetic equation [T] is a linear partially differential equation that can be solved by 
the method of characteristic. The characteristics of equation [H coincide with the ballis- 
tic trajectories of individual H atoms. The distribution function /(r,w) changes along a 
characteristic as following 

df{r,yf,t) r,( \ ft \ 
= -^(^)-/(r,w). 

After integration of the last equations taking into account the outer boundary conditions 
one can write the solution of equation [Has follows: 

/(r,w) =/5(r,,w,)-exp|^-y" fi{r)dt 

= /5(r,,w,)-exp(^- j^dt^ (2) 

Here fsi^s^'^s) is the velocity distribution function of the H atoms at the outer boundary; 
r^jWy are the radius- vector and velocity vector at the intersection of the corresponding 
trajectory (characteristic) with the outer boundary sphere. The integration is taken along 
the trajectory of atom (i.e. along the characteristic). The hot model assumes the boundary 
conditions in form of the Maxwellian velocity distribution that can be written as: 

In the hot model case the number density (n/f ), bulk velocity (V) and temperature (r^v) 
are constant at the boundary sphere. In the calculations presented here the parameters 
were taken from the results of Baranov-Malama model at the point of crossing the 90 
AU sphere and the z-axis (i.e. at the 90 AU sphere in the upwind direction). Described 
above approach corresponds to the approaches of the models that were employed for the 
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analysis of the backscattered Ly-alpha data. Such an approach takes into account of the 
heliospheric interface in the zero order of approximation. As it will be shown below such 
approach is not enough to get accurate results. 

Comparison of the hot model and Baranov-Malama model 
results 

The comparison of the results obtained in the frame of the described above hot model with 
the results of the full self-consistent two component model by Baranov and Malama (1993) 
is performed in this section in order to estimate importance of the changes of the velocity 
distribution function in the heliospheric interface on the distribution of the H atoms in the 
region inside the termination shock. Fig. [2] shows parameters of primary (populations 4) 
and secondary (population 3) interstellar atoms as a function of the the heliocentric dis- 
tance. The results are shown in the upwind and crosswind directions. As it can be seen 
from Fig. [2] that the Baranov-Malama model gives the significant differences between Tr 
temperature components at 90 AU. This is the signature of the non-maxwellian behav- 
ior of the distribution function after passing the heliospheric interface. Since the standard 
hot model operates with the Maxwellian distribution at the boundary and therefore with 
the Tr Tz components equal each to other, the Figure [21 shows 5 — 10% differences in Tr 
and components obtained in the frame of the hot model and Baranov-Malama model 
for primary interstellar atoms. This temperature component difference remains from 90 
AU upto the small heliocentric distances. For the secondary interstellar atoms (popula- 
tion 3) the distribution function is more different from Maxwellian as compared with the 
primary component. This leads to 10 — 20% difference in the number densities of the hot 
and Baranov-Malama model for this component. Differences in the temperatures for this 
population are about 15%. Note, that there are qualitative difference in behavior of 
component. This fact is connected with nonzero third moments of distribution function of 
population 3 at the boundary sphere, which does not take into account in the hot model. 

Presented plots demonstrate that there are significant differences between the results 
of the standard hot model and Baranov-Malama model. This fact leads to the conclusion 
that the heliospheric interface considerably affects the distribution of the hydrogen atoms. 
The velocity distribution function and its moments in the supersonic solar wind are key 
factors to determine the spectrum and spectral properties of the backscattered Ly-alpha 
radiation inside the heliosphere. Therefore, application of the standard hot model for com- 
putation of spectrum of backscattered Ly-alpha radiation can lead to mistakes in analyzes 
of experimental data. 

Quemerais and Izmodenov (2002) have calculated the spectrum of backscattered Ly- 
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alpha radiation by employing the results of the Baranov-Malama model. Three popula- 
tions (populations 2, 3 and 4) of H atoms where taken into account in this model. It was be 
shown in the paper that the secondary interstellar atoms (population 3), which are absent in 
the classical hot model, have significant influence on the Ly-alpha spectrum. Existence of 
populations 2 and 3 results in the asymmetry of the spectra with respect of its maximum. 
This emphasizes necessity to take into account multi-component non-Maxwellian nature 
of the velocity distribution function of the H atoms. In addition to that for more precise 
analysis of the experimental data it is necessary to employ more detailed model for the 
H atoms distribution. This model should take into account non-stationary and 3D effects 
connected with dependence solar wind parameters on time and heliolatitude. The time- 
dependent and heliolatitudinal effects were considered by Bzowski et al. (1995, 2008) 
and Bzowski, Rucinski (1995). However, these models assumed Maxwellian distribution 
function at the outer boundary that is inaccurate. In the next section new model of the in- 
terstellar H atom transport in the supersonic solar wind is presented. This model combines 
the classical hot model approach with the results of the self-consistent Baranov-Malama 
model. Therefore, in the frame of this model it becomes possible to take into account 
both the effects of the heliospheric interface and the effects connected with temporal and 
latitudinal variations of the solar wind. 



The advanced hot model of the hydrogen atoms motion into 
the hehospheric termination shock 

In order to take into account the effects of the heliospheric interface in the framework of 
the hot model, we use especial boundary condition for the velocity distribution function 
at outer sphere of 90 AU. This boundary condition takes into account Baranov-Malama 
model results. By doing this we simulate distribution functions of different populations of 
H atoms after they passed through the heliospheric interface region. 

From Baranov-Malama model we obtained values of zero, first and second moments 
of the distribution function. The moments are the number density nn, bulk velocity Vr, V^, 
kinetic temperatures Tav, Tr, and the correlation coefficient y,^; at outer boundary. Note, 
that these parameters depend on angle 6 that is angle counted from the z— axis. They are 
also different for different populations of the H atoms. Correlation coefficient is defined 
as following: 

y^(r) = -^- [ f{r,w){VR-WR){V,-w,)dw. 
nH[r) J 

Here w is the velocity of individual atom. Dependence of the number density, velocity 
and temperatures of the H atom population on the angle 6 is demonstrated in Fig. [31 Note, 
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Figure 3: Moments of the velocity distribution function calculated in the frame of the 
Baranov-Malama model as functions of the angle (0) that is measured from the upwind 
direction. All parameters except temperature are presented in dimensionless units. X- 
axis in these plots is the angle 6 in degrees. Solid curves in plots A-D correspond to 
the primary interstellar atoms (population 4), dashed curves correspond to the secondary 
interstellar atoms of population 3. Plots E and F present temperatures: solid line - Tav, 
dash-dotted line - Tj^ and dashed line - Ty for the primary and secondary interstellar atoms, 
respectively. 



that by definition: 3Tav = Tr + T^ + T^), where Tr,!^, T(p are the temperature components 
in the cylindrical coordinate system. From the axial symmetry of the considered problem 
we have V(p = and as well as correlation coefficients Vr(p = V^cp = 0. In this work to get 
the H atom parameters at the outer boundary we run the Baranov-Malama model with the 
following boundary conditions in LISM: number densities of protons and neutral hydrogen 
atoms are rip uc = 0. 06cm^^ .nn^uc = 0.18 cm^-^, respectively, relative velocity of the 
LISM is Vijc = 26.4 /, and temperature is Tuc = 6519 ^. For the solar wind parameters 
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at the Earth's orbit we assume the following numbers: Up^s = 6 cm ^, = 441.9 km/s, 
and Mach number is assumed as Me = 4.034. 

Analytical formulas for boundary condition 

Different analytical formulas for the velocity distribution function with the given moments 
can be employed. In our calculations we used three different cases: 



1. Local-Maxwellian distribution function: 



with nH,y,Tav depending on the angle 0. 

2. Three-temperature Maxwellian distribution with analytical formula that takes into 
account three kinetic temperature components (TR,T(p,Tz) 



z 



fs{r,yf) = ^-exp - 2 + 2 + 2 



l2kTR 2kT^ IkT^ 

CR=\I = W W 

V y ruH V ruH 

3. Analytical formula for three-dimensional normal distribution that employs all zero, 
first and second moments of the distribution function including the correlation coef- 
ficient (Yzr)- In this case the expression for the distribution function can be written 
as following: 

/. (r, w) = . ■ exp (- \ ( {Vr -WRf + 



-2 



5^(V«-.«)(V,-wj)) 



Here Dr = —TR,D(n = —T(p,D^ = -^T^ and k is the Boltzmann constant, thh is 
the mass of one hydrogen atom. 
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In the latter case the outer boundary condition takes into account all second moments of the 
velocity distribution function, i.e. this case provides the most complete information on the 
velocity distribution function in comparison with other cases. In the rest of this section we 
will call model 1, model 2 and model 3 is accordance with the boundary condition case 1, 
2 and 3, respectively. Fig. |4] demonstrates the comparison of the results of the models 1, 2 
and 3 with Baranov-Malama model results. As it can be seen from the plots, results of the 
simplest (first) model that takes into account dependence parameters at boundary sphere on 
angle 0, but does not take into account non-Maxwellian nature of the distribution function, 
had the worst (among the three models) agreement with Baranov-Malama model's results. 
This fact emphasizes again that the distribution function of the primary and secondary in- 
terstellar atoms are substantially different from the Maxwellan even as far from the Sun as 
the termination shock, and these differences are noticeable in the distribution of atoms at 
small heliocentric distances. Note, that the results of model 2 and model 3 are practically 
coincide each with other. Therefore, results of model 2 are shown only for the distribution 
of the velocity for the secondary interstellar atom population. These results are shown 
because for there are about 5% difference between results of the models 2 and 3. This 
difference can be explained by taking into account nonzero correlation coefficient in the 
third model, which especially affects the bulk velocity value. As it can be seen from all 
presented plots, there is a pretty good agreement between the results of Baranov-Malama 
model and the model 3. However, despite on the good agreement for all other parameters 
the behavior of the temperature (for population 3) with the heliocentric distance in the 
upwind direction is qualitatively different in the case of Baranov-Malama model as com- 
pared with all hot type models 1, 2, and 3. This means, that even knowledge of all second 
moments of the distribution function is not enough to get good adequate description for 
the non-Maxwellian nature of the distribution function for the secondary interstellar atoms 
of population 3. We suggest that this qualitative difference in the behavior of is mainly 
connected with the asymmetry of the distribution function at the outer boundary due to 
the selection effects that was discussed in Introduction. The realistic velocity distribution 
function of H atoms in the vicinity of the TS has nonzero third moments, but this fact 
is not taken into account in the considered above cases. In the next section another step 
toward the most complete and realistic description of the velocity distribution function at 
outer boundary will be done. 

On the difficulties by using the distribution function from Monte- Carlo 
calculations as the outer boundary conditions of the hot-type models 

For the most correct determination of the distribution function at the TS that takes into 
account all effects of the heliospheric interface we can in principle use the velocity dis- 
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tribution function calculated directly in the frame of the Baranov-Malama model. Indeed, 
the Monte-Carlo method with splitting of trajectories (see Malama 1991) that is employed 
for the Baranov-Malama model enables to calculate the velocity distribution function as 
it was demonstrated by Izmodenov et al. (2001). However, in order to get the velocity 
distribution with high order of accuracy that is necessary for the boundary condition of the 
hot model requires a very large number of the the trajectories in the Monte-Carlo meth- 
ods and it requires very expensive (in time and power) computer simulations. Also the 
limitations of computer operation memory becomes a strong restriction. That is why the 
Baranov-Malama model can not be directly used for the calculation of the detailed veloc- 
ity distribution function in the entire heliosphere. Nevertheless, for the calculations of the 
backscattered Ly-alpha radiation spectral properties it is necessary to know the distribution 
function of H atoms. The knowledge of the moments is not enough. 

In addition to the models 1-3 discussed above we will consider model 4 that employs 
the velocity distribution function calculated in the frame of the Baranov-Malama model by 
Monte-Carlo method. In this case the distribution function at outer boundary is obtained 
numerically and is provide in the nodes of mesh, i.e. we know values of the velocity distri- 
bution function for the given values of angle 6 and for the given three components of the 
velocity Wr,WQ,W(p. The components are in the spherical coordinate system. Using of a 
grid of uniform distribution of the grid nodes in the velocity space leads to the insufficient 
results, because only the atoms with small angular component of velocity can reach vicin- 
ity of the Sun. Therefore, small numerical errors in the interpolation in the velocity space 
at the outer boundary may lead to the essential errors near to the Sun. This explains why 
we choose a nonuniform mesh in the velocity space. This grid has smaller cells when the 
angular components of the atom velocity are close to zero. Results obtained in the frame 
of the model 4 are presented on Fig. HI It can be seen that for model 4 the T^ behaves 
qualitatively similar with the Baranov-Malama model. Let us remind that for the models 
1-3 we had qualitative difference in the distribution for the secondary interstellar atoms 
as compared to the Baranov-Malama model. Despite that behaves better for model 4, 
all results obtained in the frame of model 4 are somewhat different for Baranov-Malama 
results. These are due to the quantitative errors connected with insufficient accuracy in the 
calculations of the distribution function by Monte-Carlo method, and essentially with the 
inaccuracies in the interpolation of the distribution function at the outer boundary. Note 
that for the primary interstellar atoms (population 4) the results of model 3 have better 
agreement with Baranov-Malama model than for the model 4. This is connected with the 
fact that the distribution function of the primary interstellar atoms far from the Sun is quite 
close to the three-temperature Maxwellian distribution function. Also model 3 with an- 
alytical boundary conditions has no problems with interpolation errors, which appears in 
model 4 with tabular distribution function at the outer boundary. 
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In conclusion of this section we summarize that: 1) use of the third advanced hot 
model (model 3) with analytical formula ([3]) for the boundary conditions, which takes 
into account all second moments of distribution function, leads to acceptable agreement 
with the Baranov-Malama model; 2) for more precise accounting the global effects of 
the heliospheric interface it is possible to use realistic distribution function calculated by 
Monte-Carlo in the frame of Baranov-Malama model. However, in this case there are 
problems with accuracy of calculations. The problems are connected with the interpolation 
of the distribution function at the outer boundary between mesh's nodes in the angle 6 and 
in the three components of velocity and also due to statistical errors of the Monte-Carlo 
method. All these facts do not allow to achieve good agreement of model 4 results with 
the results of the Baranov-Malama model despite the pretty good qualitative agreement of 
the model results. 

Conclusions 

In this work we presented the comparison of the results obtained in the frame of the stan- 
dard hot model that does not take into account effects of the heliospheric interface with 
the results of the kinetic-gasdynamic Baranov-Malama model. Some possibilities for ad- 
vancement of the hot model are presented and discussed. In order to improve the hot-type 
model we need to employ the specific boundary conditions at the outer boundary far from 
the Sun. These boundary conditions should take into account the changes of the distribu- 
tion function of hydrogen atoms in the heliospheric interface region. Different cases of the 
analytical boundary condition were considered. Our main conclusions are following: 

1. The modifications of the velocity distribution function of the H atoms in the helio- 
spheric interface region has essential influence on the distribution of the interstellar 
H atoms in the heliosphere. These heliospheric interface effects leads to the sub- 
stantial differences between the results of the hot and kinetic-gasdynamic Baranov- 
Malama models. 

2. It was shown in the paper that the advancing the hot model by using local Maxwellian 
distribution function at the termination shock (i.e. Maxwellian distribution with 
parameters, which depends on angle 0) is not enough for correct inclusion of the 
heliospheric interface effects. 

3. It was shown that model 3 with the outer boundary condition ([3]), which takes into 
account all zero, first and second moments of the velocity distribution at the outer 
boundary, provides the best agreement with the results of the Baranov-Malama 
model. Main difference between the results of the model 3 and Baranov-Malama 
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model is in the qualitative difference in the temperature component of the sec- 
ondary interstellar atoms in the upwind direction. This means than knowledge of 

all zero, first and secondary moments of the distribution function at the termination 
shock is not enough to obtain complete adequate description of non-maxwellian 
character of the distribution function of the secondary interstellar atoms of popula- 
tion 3. 

4. Using velocity distribution function obtained numerically in the framework of Baranov- 
Malama model does not allow to get good agreement because of computational re- 
strictions connected with the precisions of the Monte-Carlo method and accuracy 
of the interpolation procedure. However, unlike the models 1-3 the model 4 pro- 
vides qualitative agreement of component distribution for the secondary interstel- 
lar atoms with the results of the Baranov-Malama model. 
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Figure 4: Comparison of the results of Baranov-Malama model with the advanced hot 
models. A: secondary interstellar atoms (population 3); B: primary interstellar atoms (pop- 
ulation 4). Solid curves represent the results of the Baranov-Malama model, dashed-dotted 
curves show results for the model 1, dotted curves for the model 3. Small dotted curve in 
first plot in panel A shows results for model 2. Dots in panel A correspond to the results 
of model 4. Description of different models is given in the text. 
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